arXiv:1505.04572vl [cond-mat.stat-mech] 18 May 2015 


Large deviation function of a tracer position in 
single file diffusion 

Tridib Sadhu 1 and Bernard Derrida 2 

1 Philippe Meyer Institute for Theoretical Physics, Physics Department, Ecole 
Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex 05 - France 

E-mail: tridib.sadhu@lpt.ens.fr 

2 Laboratoire de Physique Statistique, Ecole Normale Superieure, Universite Pierre 
et Marie Curie, Universite Paris Diderot, CNRS, 24 rue Lhomond, 75231 Paris Cedex 
05 - France 

E-mail: derrida@lps.ens.fr 

Abstract. Diffusion of impenetrable particles in a crowded one-dimensional channel 
is referred as the single file diffusion. The particles do not pass each other and 
the displacement of each individual particle is sub-diffusive. We analyse a simple 
realization of this single file diffusion problem where one dimensional Brownian point 
particles interact only by hard-core repulsion. We show that the large deviation 
function which characterizes the displacement of a tracer at large time can be computed 
via a mapping to a problem of non-interacting Brownian particles. We confirm 
recently obtained results of the one time distribution of the displacement and show 
how to extend them to the multi-time correlations. The probability distribution of 
the tracer position depends on whether we take annealed or quenched averages. In 
the quenched case we notice an exact relation between the distribution of the tracer 
and the distribution of the current. This relation is in fact much more general and 
would be valid for arbitrary single file diffusion. It allows in particular to get the full 
statistics of the tracer position for the symmetric simple exclusion process (SSEP) at 
density 1/2 in the quenched case. 
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Motion in crowded one dimensional channels appears in many physical systems. For 
example, passage of ions through a narrow pore in a cell membrane [I], sliding of proteins 
along DNA sequences |21, transport in ionic conductors j3], molecular motion inside a 
porous medium jij or inside carbon nano-tubes [5], etc. What makes the single hie 
diffusion interesting from a theoretical perspective is that the motion of individual 
particles is sub-diffusive. The variance of the displacement x(t) of a tracer particle over 
a time t scales as \/t as opposed to the linear dependence one finds in normal diffusion. 
The tracer particle is indistinguishable from the other particles: it only carries a mark. 
This sub-diffusive scaling has been demonstrated in several experiments on single hie 
systems (transport of water molecules inside carbon nano-tubes [B], colloidal suspensions 
in a narrow groove |&|-8], large molecules in a porous medium like Zeolite [9], etc.). 

There is a long history of theoretical works on single hie diffusion. One of the earliest 
known results is by Harris 1101 where an explicit formula for the variance was reported 


for a system of impenetrable Brownian point particles. Later, Levitt 11 generalized 
the work for a system of hard rods and derived an analogous formula for the variance. 
An interesting feature was noted in 10 -13 that at large times the huctuations of the 
tracer position become Gaussian. This was confirmed from an exact solution of the 
probability distribution of the tracer position for Brownian point particles, which is 
valid even at finite times 14 . The analysis for a system of point particles is simpler 


compared to a more general single hie problem because in the former system the 
trajectories of the impenetrable particles can be mapped to the trajectories of non¬ 
interacting particles [10|. For more general inter-particle interactions, less results are 
available: for symmetric simple exclusion process, the variance and the fourth cumulant 
are known 13,15 16 , while for colloidal systems an expression for the variance was found 
in terms of isothermal compressibility 17 . An expression of the variance for general 


single hie system was derived in 15,16 using the macroscopic fluctuation theory 18-21 


There is an interesting connection to interface huctuations where the displacement of a 


tracer can be mapped to the height of a one-dimensional interface 22 


In recent years, attention was drawn to an intriguing property of the single hie 
diffusion where the variance depends on the choice of initial state, even at large times 
There are two types of initial state one can consider 


15,16,23,24 


25 : annealed and 


quenched. In the annealed case, the starting arrangement of particles are drawn from an 
equilibrium state of average density p. This includes conhgurations where the density 
prohle fluctuates significantly from the average value. In the quenched case, only those 
initial conhgurations are considered which at a macroscopic scale have a uniform density 
prohle p (for example, a conhguration where particles are placed at uniform separation 
p ~ l ). It was found 23 that the variance is different in the two settings, although the 
y/t scaling is the same. This is surprising as one would naively expect that a system 
forgets its initial state after a long time. 

Drawing analogy with disordered systems 25 , the two settings can be described in 
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terms of cumulant generating function. For the annealed setting the generating function 
of the cumulants of the position x(t) of the tagged particle can be defined as 


M a (X,t) = In 



evolution+initial- 


(la) 


where A is a fugacity parameter and the average is over both the stochastic evolution 
and the initial state. On the other hand, for the quenched setting, the same generating 
function reads 




In 



evolution- 


^ initial 


(lb) 


where the average inside logarithm is over the stochastic evolution, whereas the average 
outside is over the initial state. Cumulants are obtained as usual by expanding M a ( A, t ) 
in powers of A: 



d n M a (X,t ) 


d\ n 


A=0 


( 2 ) 


where we use subscript a to denote the annealed and the quenched case. This notation 
will be used in the rest of the paper. 

It was found in [l5] that not only the variance, but all the cumulants of the tracer 
position scale as \fi , at large time. This implies that the cumulant generating functions 
M a (X,t ) and M q (X,t ) grow as \/t. This means that for x(t ) ~ \Jt , the probability 
distribution P a ( x(t )) of the tracer position x(t) at time t is of the form 


P a (tfVt) x ^ (3) 

where 4> a (0 is the large deviation function. The symbol Ax B denotes that the ratio 
of the logarithms of A and B goes to 1 at large t , he., lim^oo = 1. The large 
deviation functions are related to their corresponding cumulant generating functions by 
a Legendre transformation 


M a (X,t ) = Vt max 





( 4 ) 


Recently, the large deviation function </>„(£) was computed for a single hie system 
of Brownian point particles using the macroscopic fluctuation theory 15,16 and also 


In both approaches the analysis involves 


starting from microscopic dynamics 16,26 
rather long calculations. 

In this paper, we present an alternative derivation of the large deviation function 
which in our opinion is simpler. Our analysis is for impenetrable Brownian point 
particles and relics on a connection between the trajectories of the impenetrable particles 
and the trajectories of non-interacting particles, he., Brownian particles allowed to cross. 
Using this mapping we calculate the large deviation function in the annealed and the 
quenched setting, confirming earlier results (l5p6,26 1 . The large deviation functions are 
different in the two settings; even the asymptotics of the functions differ significantly: 
in the annealed setting (p a (0 ~ |£|, whereas in the quenched setting <f) q (i ■) l£| 3 , for 
large f. 
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An advantage of our method is that it can be easily extended to analyse the 
correlations of the tracer position at multiple times. The joint probability P(x(ti),x(t 2 )) 
of the tracer position x(t\) at time t\ and x(t 2 ) at time t 2 has a form analogous to 


Pa (t(Ti)Vi,Z(r 2 )Vi} x e-^ 


(5) 


where tj = rfi and x(tf) = We will show how to calculate the large deviation 

function </> a (£(Ti), £(t 2 )) for both quenched and annealed cases. This will allow us to 
recover the two time correlations of the tracer position [231 


(x{tf)x{t 2 )) = 


n(p) Vh + Vh ~ \Ah ~ h 

k(p) Vh + t 2 - vl t 


t\ 


for annealed 
for quenched 


( 6 ) 


where the pre-factor n(p) depends on the density p, and the variance a 2 /2 of the 
displacement of a single particle in unit time when it is diffusing freely in absence 
of the other particles 

cr 


«(p) = 


2 Pa/T‘ 

The time dependence of the correlation 


(7) 


in the annealed case has the same form as 


in a fractional Brownian motion 27 corresponding to the Hurst exponent H — 1/4. 

The two time correlation function calculated for more general single file problems 
using the macroscopic fluctuation theory [28] gives the same dependence on time as in (J6]) 
and only the pre-factor k(p) depends on the specific details of the system. Similar time 


dependence was found in a related model where the order of particles is preserved 29 
and also in height fluctuations of a one-dimensional interface 


Our method can further be extended to analyse the joint probability distribution 
of tracer positions at an arbitrary number of times. By taking a continuous limit this 
leads to a path integral formulation of the probability of a trajectory of the tracer as a 
functional P[x(t)] of the trajectory x(t). (We use a notation [] to denote a functional.) 
For an evolution in a time window [0, t] this probability has a form analogous to (J3]) and 
§): 


P n 


£(r)Vt 


3 —\A<M£(t)] 


( 8 ) 


where 0 < r < 1 is the time rescaled by t and 0 q [£(t)] is the large deviation functional. 
We shall show that the 0 a [£(r)] can be formally expressed in terms of the solution of a 
Schrodinger like equation with a moving time dependent step potential. 

A direct consequence of (J8]) is that the joint probability distribution of the tracer 
position at multiple times is asymptotically a multivariate Gaussian distribution. This 
can be seen by a Taylor expansion of (j) a [£(t)] around the minimum. This implies that 
all multi-time correlations of the tracer position can be expressed in terms two time 
correlation (|6]) using the Wick’s theorem (or the Isserlis’ theorem). 

Our analysis will be presented in the following order. In section [2] we define the 
problem and outline our method of calculating the large deviation function. In section [3] 
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Figure 1 . A sample trajectory of five Brownian point particles with hard-core 
repulsion. 


we present an explicit calculation of the large deviation function 4> a (0 and 4> q (£,). In 
section [4] we use this result to calculate the cumulant generating function and derive 
expressions of the first few cumulants. In section [5] we extend our method to the two 
time case, and use this to derive (|6]) . A further generalization to multi-time statistics 
of tracer position is presented in section [6j We conclude in section [7] summarizing our 
results. In the Appendix we derive a relation between the distributions of the position of 
the tracer particle and of the time-integrated current in the quenched case. This relation 
is valid for arbitrary single hie problem. It allows us to obtain the exact distribution of 
the displacement of tracer (through M q (X,t)) for the symmetric simple exclusion process 
at density p = 1/2 in the quenched case. 

2. The mapping of single file diffusion to non-interacting particles 

We consider a system of Brownian point particles on an infinite line diffusing with 
hard core repulsion between particles. A sample trajectory of the particles is shown in 
figure [l] 

At times much smaller than the mean collision time the tracer diffuses normally. 
With increasing time, the displacement of the tracer is significantly confined by the 
presence of other particles and the motion becomes sub-diffusive. We analyse the tracer 
displacement in this large time limit. 

Mapping to non-interacting particles: To derive an expression for the probability ([3]), 
we use a mapping to a system of non-interacting particles where trajectories can cross 
each other. For simplicity we discuss the mapping for a system with a finite number of 
particles which are ranked as {..., —2, —1, 0,1, 2,...} according to their positions. One 
can easily generalize the mapping to the case of infinitely many particles. In both the 
single hie and the non-interacting systems we consider the same initial condition, with 
the zeroth rank particle at t — 0 being at the origin. In the single hie system we choose 
the zeroth rank particle as the tracer and define P(x(t )) as the probability to hnd the 
tracer at position x(t) at time t. In the non-interacting system we re-rank the particles 
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at time t according to their positions and define P m (x(t)) as the probability that the 
particle of rank zero is at position x(t). These two probabilities are same 

P(x(t)) = P m (x(t )) (9) 

irrespective of the choice of the initial state as long as it is the same for both the systems. 

This equality is simple to argue from the fact that for every history of the non¬ 
interacting particles it is always possible to construct a valid history for the impenetrable 
particles, by simply interchanging the particle labels when two trajectories cross in the 
former case. As the particles are point particles, this relabelling does not change the 
probability of the history. 
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The large deviation function: In the limit of infinitely many particles distributed with 
uniform density p, the probability P(x(t )) has the large deviation form (J3]) , and due to 
(J9| the probability P m (x(t)) has also the same form. We now show how to calculate the 
large deviation function </>„(£) by analysing the large time asymptotics of P m (x(t)). 

To begin with the calculation we define Rfix) as the number of non-interacting 
particles which start at any position < 0 and reach a position > x at time t. Similarly, 
R[(x) is the number of non-interacting particles which start at any position > 0 and 
reach a position < x at time t. The quantities Rt(x) and R'fix) are random variables 
and depend both on the stochastic evolution and the initial condition. Given that the 
zeroth rank particle at t = 0 is at origin, we see that the cumulative probability of the 
zeroth rank particle at time t to be at a position x(t) > x is equal to the cumulative 
probability of the difference Rfix) — R! t {x) to be strictly positive. 


P a (x(t ) > x) 


Rt(x) 


R' t (x ) > 1 


for any x, 


( 10 ) 


where a denotes the quenched or the annealed initial condition. 

As the particles are diffusive it is expected that the probability density of the 
difference Rfix) — R'fix) has a large deviation form 


P, 


RfifVi) - RfifVt) 

Vi 


= r 


D -Vt ll>a(r,£) 


( 11 ) 


where f> a V, £) is the associated large deviation function. Then from (10) we are going 
to show that 


MZ) = M o,0 (12) 

for positive £. For negative values of £ one can use that the 0 a (£) = 4 > a (—£) which is 
expected as the microscopic dynamics is unbiased. The relation ( [l2| ) applies for both 
settings of initial condition: annealed and quenched. 


One way to derive the equality (12) is by taking the limit 

dz 


lim 

£->-oo 


in/' 


(x(t) >tvg = Hm 

£^-oo 


Vi 


Vi 


= -0a(O 


(13) 
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for positive £ where in the last step we have assumed that for positive z the 4> a {z) is 
minimum at z = 0 and is monotonically increasing away from the minimum. Similarly 
we can show that 


lim 

£—> oo 


InPjRt^Vi) - > 0 


Vi 


= -1pa( 0,0 


(14) 


where we assumed that if> a (r, £) given in ( [IT| has minimum at a negative value of r when 
£ is positive and the function is monotonically increasing away from the minimum. This 
is consistent with the property that at a position x > 0, it is more probable that Rfix ) 
is less than R[(x). Combination of the above two equations (13) and (14) along with 
(10) leads to the result (12) for any positive £. 

Rather than the large deviation function fi a (r, £) it is comparatively easier to 
calculate its Legendre transformation which gives the cumulant generating function of 
the difference Rfix) — R! t (x). Following the discussion in the introduction, the cumulant 
generating function can be defined in two settings. In the annealed setting, it is 


X JB ,?) = lim 

i-tooy/t \ / e 


(15 a) 

/ evolution+initial 

where B is a fugacity parameter and the average is over all initial configurations of 
particles and their stochastic evolution. On the other hand in the quenched setting, it 
is defined as 


£ t \ ' ' evolution j . 


(156) 


In both settings, the cumulant generating function is related to the corresponding 
large deviation function by a Legendre transformation 

Xa(B, £) = ma x{B r - ^ a (r, £)}. (16) 

r 

This leads to a parametric solution of if) a (r,£) 

dXa(B,0 


^ Q (r,£) = Br- xV B ,0 


with 


r = 


dB 


Setting r = 0 and using the relation (12) we get 

0a(£) = ~Xa(B,0 where B is determined from 


dXa(B,Q 

dB 


= 0. 


(17) 


(18) 


So our approach consists of calculating Xa{B,fi) defined in (15a, 156) and then to 
use (18) to determine the large deviation function 0 Q (£) of the position of the tracer 
particle. 


3. Calculation of the large deviation function 0 a (£) 

We now calculate Xa,{B,£) and y g (R,£). A key remark is that the difference Rfix) — 
R[{x) remains unchanged when one maps the single hie problem to the non-interacting 
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Brownian particles. For non-interacting particles located initially at positions yj, we 
can write using the independence of the particles 


,B(R t (x)-R' t (x))\ 


' evolution 


= / e BR t (x) 


-BR' t (: e ))\ 


/ evolution 


Y[ F t{Vj,B,x) (19) 


where the product is over all particles and F t (y,B,x ) is the contribution of a single 
Brownian particle starting at y given by 

d B Q(z(t)—x )' 


Ft(y,B,x) = 


( 20 ) 


for y < 0 

z(0)=y 

,- B e(x-z(t»\ for y > 0 

/ z(g)=y 

where @(x) is the Heaviside step function and z(t) is the position of the Brownian 
particle at time t. 

In this paper we consider that the initial positions of the particles are distributed 
according to a uniform density p. The annealed and the quenched cases differ by how 


the average over the initial positions is taken in the generating function (15a) and (156) 


In the annealed case, the average over initial positions is inside the logarithm. We write 
the average 


0 B(R t {x)-R’ t {xj) 


evolution+initial 


H (l -pdyj +pdyF,(y, B,x) 


( 21 ) 


where each infinitesimal interval dy is occupied by a single particle with probability p dy 
or empty with probability 1 — p dy. Taking logarithm on both sides leads to 


D B(R t (x)-R’ t (x))\ 


( 22 ) 


In ( \ = p dy \F t (y , B, x) — 1 . 

\ / evolution+initial ^ OO 

On the other hand, in the quenched case, the average over initial positions is outside 
logarithm. We get 


In (e 


D B(R t (x)-R' t (x)) 


= p dy In F t (y,B,x). 


Substituting this in the formula (15a) and (156) we get 


Xa(B,£) = p lim / dp F t (pVt,B,^Vt) - 1 
t^oo J _ 00 L V J 

XqiB, 0 = p Urn f dp In F t ( pVt , B, 


(23) 

(24 a) 
(246) 


where we defined p = yt ~2 and ^ = xt~ 2 . 

For an explicit calculation of the function F t (y,B,x ) we define a diffusion 
propagator for a Brownian particle starting at y: 


9(z,t\y,0) = 


a J na 2 t 


exp 


{z - yf 

a 2 t 


(25) 
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where ct 2 /2 is the variance of the displacement of the particle in unit time. Using this 
g(z, t\y, 0) we calculate the averages in ( [20] ) and find that F t (y, B, x ) has a scaling form 

for r/ < 0 

f(v, ~B, -0 for y > 0 
where the scaling function f(y,B,£) is given by 


F t (r)Vt,B,£,Vi ) = 


f(v,B,£) = 1 + 

with the erfc(x) being the complementary error function 


(yji.„.(s±3) 


2 f c 

erfc (a;) = / 

Jx 


du e 


Substituting this in (24a, 24 b) we get 


Xa{B,£)=pJ^ dy{f(y,B,£) - l} +p dp ^f(y,,^B, -f) - l} 

poo poo 

X,(B,()=P 'I'l I" .fin, VO + P drj 


(26) 


(27) 

(28) 

(29a) 

(296) 


Large deviation function 


Annealed: Using the expression (29a) into (18) leads to a parametric formula of the 
large deviation function 


-i) 


MO = -P—Si d ' ledc (a)~ P 

where the parameter B is determined from the relation 

9M0 




dye rfc^—j (30 a) 


dB 


= 0. 


(306) 


In this case, the B dependence can be eliminated from the formula. To show this 
we explicitly write the equation (306) using the expression of hr (30a) which leads 
to 


e B = 




dye rfcQ) ' dye rfc0) 


- 1/2 


Substituting this in the expression (30a) we get an explicit formula 

2 

r y' 


MO = \ 


J jf dp erfc 0) - dr) erfc 


a 


(31) 


(32) 


The same result was obtained previously using the macroscopic fluctuation theory 


15,16 , and also by solving microscopic dynamics 16,26 
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Quenched: Using the expression (|296[) into (18) we get 


MO = ~P dr l ln 

Jo 


1 + (f!M) erfc (i±W [ 1+ (fTW erfc (ixi) 


(33a) 


with the parameter B determined in terms of £ by solving 

dfiq(0 


dB 


= 0. 


(336) 


Unlike the annealed case, there is no explicit formula for 0 9 (£). We write the 
solution in an alternative form which will be useful later. For this we first rewrite the 


expression of <p q (f) in (33a) as 


MO =P / d V In 

Jo 


1+ (e B — l) |erfc(g) 

1 + ( e~ B — 1) |erfc(^ 3 ) 


—p / dr) ln 

Jo 


1 + 


-“*(;) 




1 + 


- erfc (^) 


We simplify the expression using an identity erfc(x) + erfc(— x) = 2 which leads to 


MO = P B M P° j d V l n 

Jo 


1 + sinlU ( — ) erfc(r/) erfc(— rf) 


where we made a change of variable rj —> a r/. With B given by (336) yields 


dri]n 


1 + sinlU ( — ) erfc(r/) erfc(— rf) 


(34 a) 


(346) 


The same formula was obtained previously using the macroscopic fluctuation theory 


15 16 , and solving microscopic dynamics 16 


Comparison of the two cases 

We see that the large deviation function in both settings of initial averaging is non- 
Gaussian. Moreover, at large values of £, the large deviation function has different 
asymptotics in the two settings, reflecting the sensitivity to the initial state, even at 
large times. To derive these asymptotics we define 


h(0 = J dr) erfc(^). 


In terms of this function the formula (32) becomes 


MO = p 2 [V h (C/o-) - VHMH 


For large positive £, the h(0 can be expanded as 

1 


h(0 = exp (—£ 2 ) 




+ 0 (£ 


-A\ 


and h(~0 = 2£ + h(0- 


Substituting this in the explicit formula (36) we get 


MO - p |4| 


(35) 


(36) 


(37) 


for large £. 


(38) 
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For the quenched case, the formula of 0 9 (£) is in an implicit form (34a, 346). Using 


the formula (34 b) we see that large positive £ corresponds to large positive B. At large 


B the integral in (34a) grows as |i? 3//2 . Substituting this in the parametric formula, 


leads to 


P 


MO * 3^2 lei 


for large £ 


(39) 


where for negative values of £ we used the symmetry £) = £> g (£). Both asymptotics 


were reported earlier in 15,16 using a macroscopic approach. 


The asymptotic behaviours (38) and (39) can be explained by a simple picture 


For the annealed case, the probability of a very large displacement x(t) of the tracer 
is dominated by the contribution from initial configurations where the interval (0, x(t)) 
is empty. The cost for the tagged particle to move a distance x(t) is then negligible 
compared to the cost of having the whole range (0, x(t)) empty in the initial condition. 
The probability of such an initial state is 


V = lim (1 — p dx 

cfcc—>0 ' 


xft)_ 

dx = e ~P*(t)' 


(40) 


Defining £ = x(t)t 2 and comparing with M we get (38). 


For the quenched case, where the initial density profile in uniform, the probability 
of large x(t) is dominated by the events where all particles inside the interval (0, x(t)) 
have reached x(t) at time t. The probability of such events, for large x(t), is 

x (t)P / /,\3 \ 

U ~ n 9 (x(t),t\p~ l j, 0 ) ~ exp \- P ^ 2 t ) ( 41 ) 


j =0 


where g(x(t),t\y,0) is the diffusion propagator (25). Using £ = x(t)t 2 and comparing 
with (J3]) we get (39). 


4. Cumulants of the position of the tracer 

The cumulant generating function is related to the large deviation function by a Legendre 
transformation Q. Using the results (34a,346) and (36) one can derive expressions for 
the cumulant generating functions in both the quenched and the annealed settings. 


Annealed 


Taking the Legendre transformation Q with the expression of £> a (£) in (36) yields a 
parametric formula of the cumulant generating function 


m „( a , t) = vif (^/m - vh-o) 

with £ determined in terms of A by solving 

x = \ (VH-O ~ VHO) erf c(£) 


erfc(£) + erfc(—£) 


h(0 


h(0 


+ erfc(—£) 


HO_ 

K-0 


h(~0 


(42 a) 


(42 b) 
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where the h(£) is defined in (35). In deriving the above formula we made a change of 

variables £ —S > <r£ in (36) and used identities h (—£) — h(£) = 2£ and erfc(£)+erfc(—£) = 2. 

The formula is simpler but equivalent to the parametric formula reported in (I5|[l6,26|. 

The cumulants of the tracer position x(t) can be extracted using (|2]). For this we 

first expand the expression of M a (X,t ) in (42a) in powers of £ and then use the second 

formula (426) to expand £ in powers of A. This can be implemented in Mathematica. 

The first non-zero cumulants obtained this way are 

, 2 \ a r 0.56419 r 

x(tf ) = ——Vt ~- ay/t 


pV 


xitfi 

x(t) 6 


o 

P 3 

a 


4 

- - 1 
7T 


P 

—j=Vt ~ 

'7T 


0.46247 


a 


P a 


Vt 


p° |_7T 


68 30 ' 

^-- + 3 


7T 


15 r 

—j=Vt ~ 

7T 


2.88197 


cr 




Vt. 


(43 a) 
(436) 
(43c) 


Quenched 

Using Q one can show that the functions 0 9 (£) and M q (X,t ) are related by 

M,(A,t) 


0,(0 = A £ 


V 


with 


1_ dM q (X,t) 

V 


(44) 


By simply comparing the above with the parametric expression (34a, 346) one can get 
an explicit formula for the cumulant generating function 


M q (X,t ) = Vo p dp In 


1 + sinlU ( — ) erfc(r/) erfc(— rj) 


(45) 


This is an even function of A which is consistent with the fact that the microscopic 
dynamics is unbiased. 

Cumulants are determined from M q (X,t ) using (J2|. As the M q (X,t ) is an even 
function of A all the odd cumulants vanish. The first non-zero cumulants are 


( x V, = 

V), = 


o 


Vt 


2 P 


dp erfc (p) erfc(— p) = 


o 


V 0.39894 


pV 27T 


p 


cr 


Vt 


(46 a) 


o 


V 


<7 2 p3 


dp 


1 — -erfc(? 7 ) erfc(— p) 
15 


erfc(r 7 ) erfc(— p) ~--- oV (466) 

P 3 


dp 1 —— 1 — erfc(r/) erfc (—p) erfc(r/) erfc(— p) erfc(? 7 ) erfc(— p) 


0.00893 


P° 


o 


v. 


(46 c) 


These expressions differ from their corresponding values (43a)-(43c) in the annealed 
case. In particular the variances differ by a factor y/2. This factor \/2 between the two 
variances is in fact present in more general single hie systems 15 16]. 


Remark: The expression (45) is very similar to the cumulant generating function of 


the time-integrated current passed through the origin |25|. In fact, one can show that 
for a general single hie problem in the quenched setting there is always a simple relation 
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between these two generating functions. We discuss this relation in Appendix A There 
is no such relation in the annealed setting. 


As the generating function of the current is known 25,311 in the quenched case 


for the symmetric simple exclusion process (SSEP) at density 1/2 we get the following 
expression for the cumulant generating function of the tracer position 

M q ( A, t) = -^= J dk In (l + sinh 2 (A)e- fc2 ) at p = \ for SSEP. (47) 

More details are presented in |Appendix A| One can easily derive the variance by 


expanding upto second order in powers of A and check that the expression is consistent 


with the known 15,16 result of variance for SSEP. 


5. Statistics of tracer positions at two times 

Our method can be extended to calculate the joint probability distribution of the tracer 
position at multiple times. As a simple example we discuss the two time case where the 
joint probability has a large deviation form (|5]). 


To calculate the large deviation function </> q (£(ti),£( 72 )) we generalize (15a, 156) 
and define two time cumulant generating functions 


Xa{B, £) = lim - 4 = In /e Ej2=1 Bj 
t-S-OO yjt \ 

Xq(B,£) = hm /in/ e ^U B o[ Rt J {xj) - R '^ Xj) ] 


evolution + initial 


evolution / initial 


(48 a) 
(48 6) 


where tj = rfi , and Xj = £,{Tj)y/t. For convenience, we defined a shorthand notation 
£ = {£(ti), £(t 2 )} and B = The R t {x) and R'fix) are defined in section [2j 

Generalising the arguments of section [2] one can show that the large deviation function 


0a(£( r i)> £( r 2 )) has a parametric solution analogous to (18): 

0 a (£(Ti),£(T 2 )) = -Xa{B,£) 
with the B determined in terms of £ by solving 

dx 


dB 3 


= 0 


for j = 1, 2. 


(49 a) 


(496) 


Then the calculation of the large deviation function reduces to determine Xa{B,£) and 


use (49a,496). 


An explicit solution requires to calculate the multi-particle averages in (48a,486) 


These averages can be expressed in terms of the contribution from a single Brownian 


particle. This leads to formulas analogous to (24a, 246): 


X a(B,£) = P hm J dp F t (r]Vi,B,£y/i) - 1 
Xq(B, £) = p lim [ dp In ' F t ( pVt , B,£\/i 


£—>■00 


(50 a) 
(506) 




























Single file diffusion 
with F t (y,B,x ) given by 


14 


Ft(y,B,x ) = 


gEj=i ©( z (b) x i) 


g Flj=i z(tj)) 


z(0)=y 


z(0)=y 


for 7 / < 0 
for y > 0 


(51) 


where the average is over a Brownian particle z(t) which started at z(0) = y. We used 
a notation x = 


It is straightforward to calculate the single particle averages in (|51|) and we present 
result. Similar to (1 

F t (riVi,B.(Vi) = 


(52) 


only the final result. Similar to (26) the F t (y, B,x) has a scaling form 

for 77 <0 

/(??, - B , -£) for 77 > 0 

where the right hand side does not depend on t and the scaling function f(rj,B , £) has 
an expression 

= 1 + h B ‘- 1 h rfcf l(Tl) d’ , i + h^W ?(T2) + " 




( U2-U l) 2 _ (•Ul+i?) 2 

(T 2 (r2—n) CT 2 ri 


/■oo /-OC exp 

+ ( eBl - !) - !) - / du ± du 2 - 2 / 7 T 

71 ^5(n) d^(r 2 ) ^ V r l( T 2 - r l) 


(53) 


Substituting (52) in (50a, 50 6| ) and then using (49a 496) we get a parametric solution 
of the large deviation function: 


For annealed case: 
0a(f(7l),f(T 2 )) = -p 


/ OO /* OO 

- i} -p J dri{f(r},-B,-£) - l| 


with i?! and h? 2 determined by solving 

dfpci n r • 1 o 

9S- = 0 fOT2 = 1 > 2 - 

For quenched case: 

/•oo /»oo 

<M£( t i)>£( t 2 )) = ~P d P 111 /( 77 , B, £) -p di 7 In/ (77, -S, -£) 

Jo Jo 

with I?! and i? 2 determined from relations 

dfi q 


(54 a) 


(546) 


as. 


= 0 for j = 1, 2. 


(55 a) 


(556) 
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The Legendre transformation of </> q (£(ti), £( 72 )) contains the information on correlations 
of all order between the tracer positions at two times. Generalizing Q we define the 
Legendre transformation 

M a (Xi,Tit\ X 2 ,T 2 t) = Vt max {Aif(ri) + A 2 £(t 2 ) - </> Q (£(o), £(t 2 ))}. (56) 

£(Ji),£(t 2 ) 

Analogous to (J2]) one can calculate the multi-time cumulants of arbitrary order by using 

d k+e M a (\i, t\] A2, t 2 ) . . 

(57) 


(xitxf x{t 2 y^j = 


d\\ dX e 2 


(Ai,A 2 )=(0,0) 


for non-negative integers k and t. Using the solutions (54a 546) and (55a, 556) one can 


then calculate the two time correlation lx{tf)x(ff))a and obtain the results (J6]) . (This 
only requires to calculate the above expressions upto quadratic order in the variables £ 
and B.) 


6. Probability of a trajectory of the tracer 

An interesting generalization would be to calculate the probability weight of an entire 
trajectory of the tracer position. For the Brownian motion this probability can be 
written as exp(—A[X]) where A[X\ is a functional of the Brownian trajectory X(t) 


which has the form of the classical action of a free particle 32,331. We want to see how 


to generalize this action for the single file problem. The probability of a trajectory of 
the tracer has an exponential form as given in Q. One can try to calculate the large 
deviation functional 0 Q [£] by generalizing the method described in section [5] 


We begin with a generalization of (51) and define a functional 


Ft[y,B,x\ = 


- f* dt'B(t') e(x(t')-z(t')) 


z(0)=y 


z(0)=y 


for y < 0 
for y > 0 


(58) 


where z(t') is the trajectory of a Brownian particle in a time window [0, t]. 


This quantity F t [y,B,x] has a simple scaling with time as in (52). We define 
77 = t~^y, £(r) = t~^x(rt), and b(r) = tB{rt ) and denote 0 < r < 1 as the rescaled 
time. In terms of these variables we find a scaling functional f[r),b,£\ such that 

./' [~V- £] for 77 < 0 
for 77 > 0 . 

Generalizing the line of arguments in the two time case we find that 4> a [f,} has a 
parametric solution in terms of f[y ,&,£]: 


Ft [», B,x] = 


(59) 


Annealed : 


( t ) a[fi\ = -pj drj^f[r], b, £] 


P 


'0 


/»oo /»oo 

Quenched-. [£] = ~p dy In f[rj, b,£]- p / dp 111 /( 77 , ~b, -£] 

Jo Jo 


, -b, -?] - q (so) 
(61) 
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where in both cases the corresponding functions b{r) are determined from a functional 
derivative 

Ha [f] 


5b 


= 0 


(62) 


with a denoting annealed or quenched. 

Unlike the earlier examples, there is no closed form solution for f[rj,b,f\. We can 
nevertheless write, by using the definition (58 59), the solution as a series in powers of 
b(r), leading to 


°° r 1 rTk -1 

f[v,*>,£] = 1 + ^2 dTl ' / 

/,■= i do 


dTkbfirf) ■ ■ ■ b{r k ) W k 


where 1-14 is given by 

m= I 


du 


l • • • 


4bi) 


4(A) 


du k g(u lt r 1 |a 2 ,’2) • ■■g(u k ,T t \ - >j, 0). 


(63) 


(64) 


Substituting this solution in (60,61) one gets the same correlations as i§ at the lowest 
non-trivial order. 

In the Appendix B we take an alternative approach and show that /[■//,&,£] can be 
expressed in terms of the solution of a Schrodinger like equation with a time dependent 
moving step potential. We did not find a way of solving this equation other than a 
perturbation expansion which also leads to (63,64). 

In 28 , the problem was approached using the macroscopic fluctuation theory. It 
was found that the Legendre transformation of </> Q [£] can be expressed formally as a 
minimum action of a variational problem, but no explicit solution was derived. 


7. Summary 

We have described here a simple method to analyse the large time statistics of the tracer 
displacement in a system of Brownian point particles with hard core repulsion. Using 
this method we reproduced the known results for the probability of tracer position at 
one time and also derived new results for the joint probability distribution of tracer 
position at multiple times. In particular, we derived a parametric formula for the large 
deviation function of the tracer position at two times and used it to calculate the two 
time correlation (|6]) . The results depend on the choice of averaging over the initial state, 
specifically the annealed and the quenched setting. In a further generalization of our 
method we took a continuous limit and obtained a formal solution of the action which 
characterizes probability of an entire trajectory of the tracer in terms of the solution of 
a time dependent Schrodinger equation. 

At large times the joint probability of tracer positions at multiple times 
is asymptotically a multi-variate Gaussian distribution with a covariance matrix 
determined by the two time correlation (J6]), instead of (z(ti)z(t 2 )) = 2D min(t 1 , t 2 ) 
for the Brownian motion. It would be interesting to find, for the single hie problem, 














Single file diffusion 


17 


analogues of celebrated properties of Brownian motion, e.g. survival probability, 
distribution of maximum, etc. 

Our method is simple, but strongly relies on the connection to non-interacting 
particles and this makes it applicable to limited systems. For example, there is no 
such connection for systems where inter-particle interactions are more complicated 
than hardcore repulsion, or when defined on a lattice (for example, symmetric simple 
exclusion process). 

In the present work we noticed an important connection between the tracer statistics 
and the current statistics in the quenched case. This connection could be useful for 
general single hie problems, in particular for the symmetric simple exclusion process 
where it has been challenging to calculate the large deviation function of the tracer 
position: since the distribution of the integrated current is known at density 1/2 for the 
quenched case 25,31 , this connection gives the full statistics of the tracer position in 
the quenched case for the symmetric simple exclusion process at density 1/2. 
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Appendix A. Relation between the distribution of integrated current and 
the distribution of tracer position 

We consider an arbitrary single hie system where the tracer always starts at the origin. 
The initial state is quenched at a uniform density p. Let x(t) be the position of the 
tracer particle at time t and let Qt(y) denotes the time-integrated current at position y 
up to time t. We define the cumulant generating functions in the quenched setting for 
these two quantities as 

n(\) = lim -^=/In /e Xx ^\ \ (1.1a) 

t t ' ' ' evolution / initial 

z/(A,£)= lim -^=/In / e x Q t( £Vt)\ \ _ (1.16) 

t tOO yj \ \ / evolution/ initial 

In fact, if the initial state has a hat prohle at density p, one finds by translation invariance 
that there is no £ dependence in z/(A,£) and 

i/(A,£) = KA,0). (1.2) 

Proposition: For a general single hie problem we are going to show that 

/'■!'A) = V (), <j) . 


(1.3) 
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Proof: We prove this by showing that the following large deviation functions are equal. 

In P(x(t) = fVt) 


m = 


lim 


t->00 y/t 


<p(q,Q = - 


f->-oo 


inP(Qt((Vt) = Wt) 
Vt 


(1.4a) 


(1.46) 


where P(x(t )) is the probability of the tracer position x(t) at time t, whereas P(Q t (y )) 
is the probability of integrated current Qfiy) at time t at position y. The y/t scaling is 
expected for a diffusive system. 


The equality between (1.4a) and (1.46) comes simply from the fact that for the 


tracer to travel a distance x(t), all the particles which were initially between the origin 
and x{t) must cross the position x{t). For the quenched initial state with uniform density 
p this number of particles is p\x(t)\. This leads to 


PCrftj) = p(q,(:r(t)) = px(t)f 


(1.5) 


Note that for the annealed case where the initial density is fluctuating, this equality 
would not hold. 


Using the above equality and the large deviation form (1.4a, 1.46) we get 

0R) = />.£)■ 


( 1 . 6 ) 


The cumulant generating functions //(A) and v{ A, £) are related to the corresponding 
large deviation functions 0(£) and <p(q,£ ;) by a Legendre transformation 

= max{A f - 0(0} (1-7) 

^(A,0 = max{A (1.8) 


Taking the Legendre transformation and using the translation invariance (1.2) we obtain 


the identity (1.3). 


An example 


We now check the relat ion ( 
repulsion. Comparing 


calculated M q (\,t) in ( 


1.1a 


5) 


3) for a system of Brownian point particles with hard core 


and (1&) we find that //(A) = t where we have 

dere we show that the cumulant generating function u{ A, 0) 


of integrated current is given by 

poo 

v(\, 0) = p a / drj In 


1 + skill 


— ) erfc(^) erfc(— rj) 


(1.9) 


The two results are consistent with (1.3). 


To derive (1.9) we realize that for the integrated current we can simply ignore the 


hard core repulsion between particles and treat them as independent particles. Following 
a calculation similar to the one in section [3] we find 
AQd o)\ 


/ evolution 


rufe,A.o) 


( 1 . 10 ) 


J 
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where F(y, \,x) is the function (20) and j denotes all particles. Using the scaling (26) 


it is simple to show that the cumulant generating function 

POO POO 

v{\ 0) = p / drj 111 /( 77 , A, 0) +P / drj In 0). 


( 1 . 11 ) 


A straightforward algebra using (27) and the identity erfc(x) + erfc(— x) = 2 leads to 


the result (1.9). 


Cumulant generating function for SSEP 

For the symmetric simple exclusion process (SSEP) the generating function A, 0) of 


the integrated current is known 25 at density p — 1/2. The expression can be derived 
by using a relation between the cumulant generating function of integrated current in 
the quenched and the annealed settings [25] Eq. 43]: 


Z/(A, 0) — Vquenchedfif'i 0) \/'2 ^ annea ^ e ^^^ ^ ^ d 2 SSEP, (1.12) 


where the latter quantity is defined as 

1 - / ,XQ t {Q)\ 


i2(A, 0 ^annealed lilU In ( e 

t^-OO yft 


/ evolution+initial 


(1.13) 


For SSEP the is a nneaied( A, 0) was calculated using Bethe ansatz in 31 . At any uniform 
density 0 < p < 1 it has an expression 

i r° .. . ( . o /A' 


V anneaied(\ 0) = - / dk In 1 + 4p(l - p) sinh - e 


7 r 


-k- 


(1.14) 


Substituting above expression in (1.12) and using the identity (1.3) we arrive at the 


result (47). 


Appendix B. A formal solution of f[r),b, £] 

To analyse this scaled functional f[rj,b,f\ in ( [59] ) we define 
G(u, t \ t ), 0 ) = 

\ / z( 0 )=i; 

for 0 < r < 1, where the subscripts and the superscripts on () denotes the hxed initial 


( 2 . 1 ) 


and final positions of a Brownian motion z(r i), respectively. From (58) and (59) one 
can see that 


f[v,b,€] = / duG(u, 1| - r), 0) 


Considering a small increment of r in the formula (2.1) we write 


POO 

G(u, t + e|j/, 0) — / dwg(w, e|0, 0)G(u + w, t \ tj , 0)( e^ T rfT i b f r O 0 ( : ( 7 'i)-s( T i)l\ 


( 2 , 2 ) 


z(r+e)=u 


/ z(r)=u-\-w 


where g(x,t\y,0) is the diffusion propagator defined in (25). 
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For small e we expand the right hand side in powers of e and obtain 

2 

G(u, r + e\r), 0) = G(u,r\rj,0) + e|6(r)@(n - £(t))G(u, r\rj, 0) + °^d uu G{u, r\rj, 0)| 


+ ... 


where we used (25) and g(w, 010, 0) = 5(w). Taking e —> 0 we arrive at 

_2 


cr 


<9 T G - —d uu G = 6(r)0 (« - f (r)) G. 


(2.3) 


This has the form of an imaginary time Schrodinger equation with a time dependent 
step potential on the right hand side. The initial condition for G can be determined 
from its definition (2.1) which leads to 


G(u, 0 |? 7 ,0) = 5(u — rf). 


(2.4) 


Solution of (2.3) combined with relation (2.2) determines the functional f[rj,b,f\. 


It is difficult to solve (2.3). For small b(t) one can however write a series solution 


of (2.3) as 


f * T k —1 


G(u, t\tj, 0) = g{u, t\tj, 0) + ^ / dr i - - - / dr k b (n) ■ ■ ■ b ( r k) 


k =1 


du\ 


U(n) 


'£(n) 


du k g(u,T\ui,n) ■ ■ ■ g{u k -i,Tk-i\u k ,Tk) g{u k ,T k \r), 0) (2.5) 


Substituting this in (2.2) leads to the solution (63). 
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